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ABSTRACT 

In light of recent findings from the kinematic morphology-density relation, we investi¬ 
gate whether the same trends exist in the original morphology density relation, using 
the same data as Dressier. In addition to Dressler’s canonical relations, we find that 
further refinements are possible when considering the average local projected density 
of galaxies in a cluster. Firstly, the distribution of ellipticals in a cluster depends on 
the relative local density of galaxies in that cluster: equivalent rises in the elliptical 
fraction occur at higher local densities for clusters with higher average local densities. 
This is not true for the late-type fraction, where the variation with local density within 
a cluster is independent of the average local density of galaxies in that cluster, and is 
as Dressier originally found. Furthermore, the overall ratio of ellipticals to early-types 


in a cluster does not depend on the average density of galaxies in that cluster (unlike 
the ratio of lenticulars to disk systems), and is fixed at around 30%. In the paradigm 
of fast and slow rotators, we show that such an elliptical fraction in the early-type 
population is consistent with a slow rotator fraction of 15% in the early-type popula¬ 
tion, using the statistics of the ATLAS 313 survey. We also find the scatter in the overall 
ratio of ellipticals to early-types is greatest for clusters with lower average densities, 
such that clusters with the highest elliptical fractions have the lowest average local 
densities. Finally, we show that average local projected density correlates well with 
global projected density, but the latter has difficulty in accurately characterising the 
density of irregular cluster morphologies. 
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1 INTRODUCTION 

It is well known that early-type galaxies (ETGs) are more 
common in denser (cluster) environments, with a corre¬ 
sponding decrease in numbers of late-type galaxies (LTGs: 
|Oemler 1974| Davis & Geller 19761. Dressier ( 1980b[ here¬ 
after D80) convincingly demonstrated this morphology- 
density relation in a quantitative manner, showing how the 
number fraction of lenticulars (SOs) and, at higher densities 
ellipticals (Es), rises with local projected density (logE) to 
replace a diminishing LTG (or spirals and irregulars, S+I) 
fraction. With morphology described by the T-type, this re¬ 
lation was coined the T — E relation. |Whitmore et al.| ( |1993| ) 
later showed that an equivalent relation holds with scaled ra¬ 


dius (the T-R relation). Dressier et al. (1997) also extended 
the study to higher redshifts and found a lack of SOs and an 
excess of Es, suggesting very different formation epochs for 
the two morphologies. 

Recently, integral-field spec troscopy (IFS) s urveys such 
as SAURON ( Bacon et al.|2001 1 and ATLAS 3D (Cappellari 


et al.|201 l) have been able to refine the S0/E division, based 
on the internal kinematic properties of the galaxies. In the 
complete volume limited ATLAS 3D sample, two thirds of 
the Es show rapid rotation in their velocity maps, which are 
indistinguishable from the velocity maps of typical SOs. This 
led to the classification of fast rotators (FRs), a population 


consistent with being oblate axisymmetric spheroids (Cap- 
pellari et al.|2007 Emsellem et al.|2011 1. The remaining Es 


(and some SOs) showed an absence of rotation in their veloc¬ 
ity maps leading to the classification of slow rotators (SRs), 


a population which may be (mildly) triaxial (Emsellem et al. 


2011 ). 


The T — E relation was revised by |Cappellari et ah] 
12011), who used the kinematic classifications FR/SR and 
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the morphological LTG classification to develop a kinematic 
morphology-density relation (kT — E). This kT — E relation 
shows a clear decline in the LTG number fraction towards 
higher projected density mirrored by a corresponding rise in 
the number density of FRs This is qualitatively similar to the 
previous results of D80 for the SO and LTG morphological 
classes. |Cappellari et ah] ( |2011[ ) also reported that the total 
fraction of SRs was constant at 4% across all environments, 
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except for the core of the Virgo cluster, where it rose to 
20%. This behaviour is noticeably different to that of Es in 
the original T — E relation, which show a more gradual and 
sustained increase in numbers with projected density with 
the fraction of ellipticals reaching 40%. 

However, the ATLAS 30 survey mostly sampled field 
and group galaxies, and contained only a single cluster 
(Virgo) which is unrelaxed, spiral rich, and x 100 lower in 
density than the nearby Coma cluster. Conversely, D80 stud¬ 
ied 55 rich clusters and a single ‘field sample’ to construct 
the original T— E relation. Thus, the ATLAS 30 kT — E and 
original T —E relations probe very different types of environ¬ 
ment and are not directly comparable. Motivated by this and 
the over-abundance of SRs in the core of Virgo, D’Eugcnio 
|et al.| ( |2013[ ) and |Houghton et al.| ( |2013| ) studied the kT — E 
relation in two more massive, denser clusters: Abell 1689 
and Coma. These studies highlighted that, while the internal 
fraction of ETGs and LTGs may vary, the average fraction 
of SRs in the ETG population (the (SR:ETG) fraction) is 
constant at around 15%. This result was confirmed in lower 
mass and unrelaxed clusters (Scott et al.|2014 Fogarty et al. 


20141 and appears to hold across five orders of magnitude 
in projected density, from the field (ATLAS 30 data) to the 
densest cluster (Abell 1689). However, within each cluster, 
the SR:ETG ratio varies strongly with logEs: more SRs were 
found in the densest (viz. central) regions, but only at the ex¬ 
pense of the less dense (viz. outer) regions, where FRs were 
abundant. This led to the conclusion that the SRs were be¬ 
ing segregated from the rest of the population, possibly via 
dynamical friction due to the massive nature of most SRs. 
Crucially, the SR:ETG ratio peaks at progressively higher 
logE 3 for clusters where the average projected local density 
of galaxies is higher (Houghton et al. |2013[ ). Indeed, at the 
highest projected galaxy density of Virgo (containing the 
highest SR fraction in Virgo), D’Eugenio et al.| ( |2013| ) found 
no SRs in Abell 1689. 

The implications of these finding are significant. The 
different kinematic structures of FRs and SRs are believed to 
originate from different assembly mechanisms, meaning the 
kinematic classifications directly relate to different assembly 
histories. One would not expect the efficiency of these two 
mechanisms to vary in the same way. It is often suggested 
that SOs result from quenched spiral galaxies while ellipti¬ 
cals originate via mergers (Toomre fe Toomre|1972' 


fe Hernquist||1992[ |Naab fc Trujillo||2006| |Bois et ai.||2011[ 


Barnes 


Khochfar et al.|2011al; the rapid rotation in FRs is strongly 


suggestive of a disk-like origin whereas the random, often 
decoupled, motions of a SR are indicative of a complicated, 
perhaps violent, assembly from multiple components. It is 
intriguing therefore that the overall density of galaxies in 
any environment feeds these two processes equally, produc¬ 
ing a fixed ratio of FRs and SRs. Is the formation of SRs 
dependent on the availability of FRs (i.e. a dry-merger sce¬ 
nario)? Could FRs and SRs in fact be produced by a single 
mechanism (i.e. degrees of harassment or minor merging)? 

In this paper we investigate if recent results from the 
kT — E relation also hold when using just visual morpholo¬ 
gies. Specifically we aim to find out: how the E:ETG ratio 
varies with logE 3 within Virgo-like, Coma-like and Abell 
1689-like clusters; and how the (E:ETG) ratio of an entire 
cluster varies with the average local density of galaxies in 
that cluster. The structure of this paper is as follows: 1 


describes the samples from which the data was compiled, 
n describes the results, which are then discussed in )j4] 
and fj5] concludes. We assume standard WMAP7 cosmology 
throughout (Komatsu et al.[2011). 


2 SAMPLES & DATA 


The data used in the original T— E was published in |Dressler| 
(1980a I. Using that catalogue, we updated the cluster red- 
shifts using the latest values from NEEQ as well as the dis¬ 
tances to the clusters, assuming all clusters are at rest with 
respect to the Hubble flow. We further checked the cluster 
redshifts using the galaxy redshifts available in NED (mostly 
from th e 6dF survey, | Jones et al.|2004j ), the 2MASS redshift 
survey (Huchra et al. 20121, the WINGS survey (|Moretti 


et al.|2014 1 and the CAIRNS survey ( Rines et al.|2003 1. We 


found that the redshifts published by NED for the clusters 
Abell 838, Abell 1631 and DC2349-28 were inconsistent with 
the redshifts of all the galaxies: these three objects appear to 
be projections of at least two clusters along the line of sight. 
Further details are included in Appendix A; here we simply 
exclude these clusters from all subsequent analysis. |Dressler] 
|fe Shectman| ( |T988| ) also found Abell 1631 to be a projection 
of multiple clusters, as well as Abell 1736; we do not reject 
Abell 1736 from our sample as the two components are only 
separated by 3500 km/s, which at face value only equates to 
a distance error of ~15% for each component. This results 
in a catalogue of 5455 galaxies in 53 clusters^] 

We generalise the more detailed morphological classes 


available in the Dressier (1980a) catalogue into LTGs (spi¬ 


rals and irregulars, S+I), lenticulars (SO) and ellipticals (E). 


The local projected densities of the galaxies were not 


published in Dressier (1980a I so we re-calculate them. We 


define the Nth average projected density Ejv to be N times 
the reciprocal of the smallest circular area (ajv measured in 
square arc seconds, or An when measured in kpc 2 ) at the 
distance of each cluster that encloses the nucleus of the Nth 
nearest neighbour. Where galaxies are found at the edges 
of the photographic plates, we correct on for the region 
missing of the plate edge. We adopt the same correction 
for foreground/background galaxies as D80 did, but instead 


1 http://ned.ipac.caltech.edu 

2 We do not combine data for DC 0107-46 and DC 0103-47, al¬ 
though considered a single cluster in D80. 



log(E 3 /Mpc 2 ) 


Figure 1. The distribution of the average projected density of 
galaxies in each cluster, log(E 3 /Mpc 2 ). Dashed lines divide the 
sample into low-, medium- and high-density clusters. 
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log(S 3 /Mpc 2 ) 


Figure 2. The distribution of local projected galaxy density 
(logS 3 ) within low-, medium- and high-density clusters. The dis¬ 
tributions are log-normal and don’t exhibit particularly strong 
tails to low or high density. Thus the mean projected density 
of galaxies in each cluster is a fair representation of the typical 
projected density of galaxies in each cluster. 


of subtracting this correction, we incorporate it into the cal¬ 
culation of aN'- we count the number of galaxies above what 
is expected by the background while increasing the circular 
area on the sky until the required number N is attained. In 
one case, this approach fails as the number of galaxies never 
exceeds the background rate; we remove this galaxy from all 
subsequent analysis. In three cases, the number of galaxies 
above the background rate never reaches N, so we calculate 
Ejv using the peak value of the number of galaxies found 
above the background rate. For all the other 5451 galaxies, 
this approach succeeds and provides a better measure of Ejv 
than subtracting a correction, a posteriori : N refers to the 
number of galaxies above the background (not including it) 
which helps produce a Normal distribution, particularly at 
lower densities (c.f. Fig. H & §- 

We separate clusters into groups of different (projected) 
density by taking the average projected density of the galax¬ 
ies in each cluster and defining low-, medium- and high- 
density clusters as those with average projected densities 
of log(E3/Mpc _2 )^1.25, 1.25 <log(E3/Mpc _J )^1.75 and 
log(E3/Mpc _2 }>1.75. The number of clusters in each of 
these ranges is 18, 30 and 8, respectively. Fig. |T] shows the 
distribution of log(E 3 ) and these divisions. In Fig.[2]we show 
the distribution of logE 3 for galaxies from clusters in each 
group. These distributions of logE 3 are Normal-like and not 
obviously skewed towards higher or lower density for any 
density class. Hereafter, logE 3 is used to refer to the pro¬ 
jected local density of an individual galaxy, while log(Es) 
is used to refer to the mean projected local density of all 
galaxies in a cluster. 

We calculate the average E, SO and S+I fractions in 


each cluster to compare with log(E 3 ). At this point, we are 
able to correct for differential foreground/background con¬ 
tamination using the fractions found in Dressier et al. (1997) 
and the areas of the photographic plates. In one case, the 
expected number of E contaminants is greater than the ob¬ 
served number of E galaxies; we therefore set the corrected 
number of Es to zero. 

In order to better understand the role of log(S 3 ), we 
compare it to a measure of global density. We determine the 
median RA and Doc of all galaxies in each cluster and adopt 
this as the cluster centre. While we investigated luminosity 
weighted and mean centres, we found the median to be the 
most robust. We then calculated the standard deviation of 
the galaxy locations along the RA and Dec axes and took 
the geometrical mean as a fiducial measure of the size of the 
cluster (hereafter oxy)- The sizes measured this way were 
always less than 1 Mpc. Counting the number of galaxies 
within a circle of radius oxy centred on the median coor¬ 
dinate, we measured the average global density of galaxies 
E CT = N(r<CTjvy) / (naxY)- While the results of this den¬ 
sity metric are qualitatively similar to an average within a 
fixed radius aperture (e.g. 1 Mpc), an adaptive region over 
which to calculate the density takes account of the rela¬ 
tive differences in cluster sizes and morphologies; that said, 
there are three clusters for which even an adaptive circular 
aperture fails to encompass the bulk of the galaxies due to 
their highly non-spherical distribution of galaxies. This is 
discussed in more detail in © 

For morphological fractions, we calculate posterior val¬ 
ues and uncertainties as in Houghton et al. (20131; the differ¬ 
ence between measured number fraction and posterior num¬ 
ber fraction is usually small unless the fraction is close to 
zero or one, or if the total number of galaxies is small. Our 
uncertainties (and error bars) mark the 15.9 and 84.1 per¬ 
centiles of the posterior probability distribution (i.e. ±lcr); 
similarly, our reported values (symbols on plots) mark the 
median (50th percentile) of the posterior probability distri¬ 
bution. Shaded regions on plots are drawn between the 15.9 
and 84.1 percentiles (i.e. they are equivalent to ±1 <t error 
bars). 


3 RESULTS 

We now present the original T — E relation, the T — E rela¬ 
tions for clusters of different average densities, followed by 


Table 1. Polynomial coefficients for the fits made to various re¬ 
lations between morphological type and projected density (/ = 
JT c;(log£ 3 )*). Relations 1-3 are shown in Fig. |3l relations 4-6 
are shown in Fig.[5]and relation 7 is shown in Fig.[6f Uncertainties 
are only given for the (inverse variance weighted) linear fits. 


Relation 

CO 

Cl 

C2 

C3 

1: f(E)-E 

0.00579 

0.02545 

0.01570 

0.11630 

2: f(S0)-£ 

-0.00952 

0.01434 

0.08544 

0.32620 

3: f(SI)-£ 

0.00462 

-0.04104 

-0.10229 

0.55909 

4: f((E))-(£) 

0.065T0.02 

0.063T0.01 

- 

- 

5: f«S0»-<£> 

0.272T0.03 

0.116T0.02 

- 

- 

6: f«SI))-<£> 

0.601T0.03 

-0.178i0.02 

- 

- 

7: f((E:ETG))-£ 

0.242T0.03 

0.026i0.02 

- 

- 

8: f((S0:Disk})-£ 

0.295T0.03 

0.185i0.02 

- 

- 

9: log(£ 3 )-log£ CT 

-0.249 

1.148 

- 

- 
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Figure 3. Reproduction of the T — S relation of D80 using up- 
to-date cluster redshifts, projected densities to the third nearest 
neighbour and field-edge or completeness corrections. 


the (T)-(E) and (E:ETG)-(E) relations to study how the 
average morphological mix in a cluster depends on the av¬ 
erage density of galaxies in that cluster. We further show 
how average local density relates to global density via the 
log(E 3 )-logE (T relation. 

Where polynomial fits are made and shown in figures, 
the coefficients are given in Table [I] Linear least-squares fits 
were inverse-variance weighted for all except the log(E 3 )- 
logEo- relation; for this we used a Bisector estimator in the 
absence of uncertainties on either axis (Isobe et al. 1990). 
We tabulate the coefficients and uncertainties in Tableland 
show la ranges for the fits as dotted lines in Figs. M 
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3.1 The original T — E relation 

Fig. [3] shows our re-processed data in the same style as pre¬ 
sented by D80. Using the number fractions of LTGs, SOs and 
Es, we recover the same trends. We fit polynomials to the 
data (unweighted by uncertainties) to act as references for 
future comparisons (coefficients in Table [I]). Note that our 
use of logEs (rather than logEio as in D80) shifts the trends 


Figure 4. The T— E relations for different density clusters. In all 
panels, we show number fraction against local projected density 
for the low-, medium- and high-density clusters (light, medium 
and dark shading; defined in Fig. |TJ) ; the black lines are the same 
as those shown in Fig. i Top: the S+I (LTG) fraction. Middle: 
the SO fraction. Bottom: the E fraction. While the variation in 
the S+I and SO number fractions is similar for clusters of different 
average density (and to the relations in Fig. |5J, the E fraction is 
dependent on both the local projected density and the average 
projected density galaxies in the cluster. 
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to higher values of local density. Furthermore, the updating 
of the sample and the removal of three clusters results in a 
higher peak E-fraction. However, the value of this peak is de¬ 
pendent on the binning: to sample a larger number of galax¬ 
ies and obtain better statistics, the bin at highest densities 
covers a larger range in local density (2.5< logE 3 ^ 3.3). 
If we split this last bin into two, the E-fraction peaks at 
a higher value, but with considerably larger uncertainties; 
in the extreme case, we could define a bin that samples 
just galaxies in the densest environments of the high-density 
clusters, thereby recovering the peak E-fraction that we will 
see in Fig. [4] 

3.2 The T — E relation for low, medium and 
high-density clusters 

Fig. a illustrates the T — S relations for low-, medium- and 
high-density clusters. Accordingly, we introduce the follow¬ 
ing terms: the f(SI)-E relation describes the change in the 
number fraction of S+Is with logEs; the f(S0)-E relation de¬ 
scribes the change in the number fraction of SOs with logS 3 
and the f(E)-E relation describes the change in the number 
fraction of Es with logE 3 . 

The maximum (viz. central) degree of segregation (i.e. 
the maximum E-fraction, maximum SO fraction and mini¬ 
mum S+I fraction) appears similar in low-, medium- and 
high-density clusters. This is contrary to what one might 
expect: clusters with galaxy populations that extend to high 
logEs appear to have similar morphological fractions in their 
densest regions as clusters whose galaxy populations only 
extend to lower logE 3 . The SO fraction also appears to drop 
at the very densest regions; this was also seen by Whitmore| 

|etUJ I pilfo 

Fig. HI also shows the best-fit polynomials from the T— E 
relation in Fig. [3] (over plotted as black lines). The variation 
of the f(SI)-E and f(S0)-E relations is broadly similar to 
these relations and all are close to linear in logE 3 . But the 
f(E)-E relation is quite different: it reaches a peak value 
of 50%-70%, irrespective of the average density of galax¬ 
ies in the clusters. Furthermore, the shape of the f(E)-E 
relation is similar for low- medium- and high-density clus¬ 
ters, but shifted to higher densities for denser clusters. Note 
that the densities which have the highest fraction of Es in 
low-density clusters (or medium-density clusters) have a sig¬ 
nificantly lower E fraction in high-density clusters. 

3.3 The (T)—(E) relation 

Fig- i shows the morphological fractions and local galaxy 
densities averaged within each cluster; individual points rep¬ 
resent different clusters while the lines show linear least- 
squared fits. The (T)-(E) relations are much weaker than 
the T — E relations, varying only slightly with the log(E 3 ) 
of each cluster. We measure the observed scatter a 0 (and 
estimate the intrinsic scatter cr,) in each relation to be 0.07 
(0.06), 0.09 (0.07) and 0.11 (0.98) for the E, SO and S+I 
fractions, respectively. The similarity between a 0 and ai 
suggests (j i is considerably greater than the measurement 
uncertainty. The scope of each linear relation (A/) over the 
range in density probed (AlogE 3 = 2.5 —0.5) is only 1.1, 1.7 
and 1.8 (for E, SO and S+Is, respectively) times the range 



Figure 5. The (T)—(£) relation, relating the average morpho¬ 
logical fraction to the average projected density (for individual 
clusters). The variation with average density is weak and compa¬ 
rable to the intrinsic scatter in each relation. 


spanned by the intrinsic scatter (±l<Ti). As such, we confirm 
that these relations are weak in their predictive power. 

Fig. © shows the mean SO fraction within the disk 
population ((S0:Disk) = N[S0]/{N[S0]+N[S+I]}) and the 
mean E fraction within the ETG population ((E:ETG) 
=N[E]/{N(E)+N(S0)}) in each cluster as a function of 
log(Es). The f((S0:Disk))-E relation shows an increase in 
the fraction of passive disk systems with log(E 3 ). This is to 
be expected given Fig.© Conversely, the f((E:ETG))-E rela¬ 
tion is consistent with a slope of zero and the (E:ETG) ratio 
appears constant at around 30%. However, there is signifi¬ 
cant scatter at lower densities ((log(E 3 /Mpc 2 )} < 1.5) with 
(E:ETG) ratios within individual clusters varying between 
10% and 50%. If we estimate the mean (E:ETG) ratio for 
low-, medium- and high-density clusters, we find 0.296+0.03, 
0.271+0.02 and 0.275+0.03, respectively (quoting the stan¬ 
dard error in the mean). We note that the standard error in 
the mean does not take account of measurement errors. 

Below log(E3/Mpc _2 )=1.5, there are a group of clusters 
with remarkably high (E:ETG) ratios. We can split these 
lower density clusters into those that have (E:ETG) ratios 
above and below 40%. The clusters in both divisions have 
similar distances (the mean redshifts are 0.051 and 0.043 for 
the high and low (E:ETG) clusters, respectively) suggesting 
the difference is not a result of a distance-dependent mis- 
classification error. Of the seven clusters with high (E:ETG) 
ratios, five have Bautz-Morgan classifications in NEE© they 
span all types (I, I/II, II/III and 2xIII). All but one of the 
clusters with low (E:ETG) ratios have Bautz-Morgan clas¬ 
sifications in NED: again they span all types (6x1, 2x1/11, 
5x11, 2xII/III and 10xIII). Thus the high and low (E:ETG) 
clusters do not appear to have substantially different mor¬ 
phologies. 

We show the T — E relations for the high and low 
(E:ETG) ratio clusters in Fig. © While the S+I fraction 
is similar in both cases, the f(E)-E and f(S0)-E relations 
appear to converge for clusters with higher (E:ETG) ratios: 
the SO fractions are lower and the E fractions higher. 


3 http://ned.ipac.caltech.edu 
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Figure 6. The (E:ETG) and (SO:Disk) ratios (red and green 
points, respectively) averaged over individual clusters, versus 
log(E 3 ). While the fraction of SOs in the disk population shows 
a steady increase with log(S 3 ), the fraction of ellipticals in the 
ETG population is constant at around 30% for all clusters. 


High <E:ETG> clusters . 

Low <E:ETG> clusters 
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Figure 7. Left: the T — S relation for low-density clusters 
with high (E:ETG) ratios. Right: the T — E relation for low- 
density clusters with low (E:ETG) ratios. All: Clusters with 
higher (E:ETG) ratios appear to have similar f(E)-E and f(S0)- 
E relations while these relations diverge for clusters with low 
(E:ETG) ratios. Note that the S+I fraction is similar in both 
cases. The black lines represent the polynomial fits in Fig. [3] 


3.4 The log(E 3 )-logE CT relation 

In Fig. [8] we show the log(E 3 )-logEo- relation, demonstrat¬ 
ing that the mean local density correlates well with the 
global density of the central cluster regions. The size of 
the points represents the fraction of all cluster galaxies con¬ 
tained within the adaptive aperture: there are three clusters 
for which the number of galaxies contained in the adaptive 
aperture is exceptionally low (<10%; Abell 548, Abell 2256, 
& DC0326-53) as the median coordinates of the cluster 
galaxies do not adequately represent the centres of the clus¬ 
ters (or substructure in the clusters). Unsurprisingly they 
are Bautz-Morgan types III, II/III and III, respectively. If 
all galaxies were distributed Normally in a circularly sym¬ 
metric distribution, we would expect 39% of galaxies to be 



logfEff/Mpc 2 ) 


Figure 8. The log(E 3 )—logS CT relation, linking the average local 
density of galaxies with the global density in an adaptive circular 
aperture. Symbol sizes represent the number of cluster galaxies 
contained in the adaptive circular aperture. The dashed line is 
the 1:1 relation, while the solid line is a bisector fit to the data. 


contained within our adaptive aperture. In practice, most 
(all but one) clusters contain fewer galaxies than this in the 
adaptive aperture, because cluster galaxies are neither Nor¬ 
mally, nor circularly symmetrically distributed on the sky. 

We fit the data using a Bisector estimator (in the ab¬ 
sence uncertainties on either axis, Isobe et al. 1990). While 
the link between average local density and global density is 
very close to the 1:1 relation (dashed line), the low-density 
clusters scatter below it, suggesting they are less well rep¬ 
resented by a circular aperture. Similarly, there is a ten¬ 
dency for global densities to be underestimated when fewer 
galaxies are contained within the adaptive aperture (at fixed 
log(E 3 ), symbols are smaller towards lower logE CT ). Further¬ 
more, clusters with higher global densities tend to have more 
galaxies in the adaptive aperture (at higher global densities, 
a larger fraction of the cluster galaxies are contained within 
the adaptive aperture). 


4 DISCUSSION 

In the following section, we argue for the following: the f(E)- 
E relation is not driven just by logE 3 , but is also depen¬ 
dent on log(E 3 ), which suggests a different mechanism is 
responsible for creating this relation; in any given cluster, 
the (E:ETG) ratio is constant at around 30% which is con¬ 
sistent with a (SR:ETG) fraction of 15%; the higher scatter 
in the (E:ETG) ratio at low average densities is real, and 
clusters with high (E:ETG) ratios all have high E-fractions; 
the log(E 3 )-logEo- relation gives meaning to our measure of 
average local density, and highlights our inability to directly 
measure an average global density. Finally we discuss the 
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implications of these findings in the context of formation 
mechanisms for ETGs. 


4.1 The f(E)—E relation is different to the f(SO)—E 
&: f(SI)—E relations 

The number fractions of S+I and SO galaxies in Fig.[f]appear 
to be independent of the log(E 3 ): we see the same f(S0)-E 
and f(SI)-E relations for low-, medium- and high-density 
clusters. This is entirely consistent with the results of D80. 

However, the situation is more complicated for the Es. 
The original T — E relation of D80 can be interpreted as a 
universal link between all morphologies and local projected 
density. Fig. [4] shows that the relation between E-fraction 
and logE 3 is not that simple: the E-fraction depends on the 
relative local density within the cluster. In low-density clus¬ 
ters, the E-fraction attains a similar peak value to the high- 
density clusters, but at appreciably lower local density. This 
peak in the E-fraction always occurs at the highest local den¬ 
sities available in a given cluster; that is, the highest relative 
local densities. 

The results from the kT — E relation are strikingly sim¬ 


ilar to the results presented here. In Houghton et al. (20131, 


the SR:ETG ratio reached peak values at very different lo¬ 
cal densities for Virgo, Coma and Abell 1689. The similar 
behaviour of the SR:ETG and E:ETG ratios suggests that 
they are tracing the same, or at least a highly correlated, 
population. This is despite 66% of the Es in the ATLAS3D 
survey being Fast Rotators, which are kinematically indis¬ 
tinguishable from SOs, and grouped with the majority of SOs 


in any kT — E relations (Emsellem et al. 2011 Houghton 
et al.|2013 1. 


4.2 The (E:ETG) ratio is around 30% in any 
given cluster 


Fig.[6]shows that, while the ratio of star-forming to quiescent 
disks (S0:Disk) strongly depends on log(S 3 ), the (E:ETG) 
ratio is indifferent to it: the (E:ETG) ratio is always around 
30%. 


This behaviour is remarkably similar to that of 
the (SR:ETG) ratio (Houghton et al. 20131. While the 
(SR:ETG) ratio was found to be to be 15%, here the 
(E:ETG) ratio is 30%. The ATLAS' 513 project kinematically 
classified 253 ETGs and we can use these statistics to link 
kinematic morphologies with visual morphologies. Let us 
define the ‘probability of error’ in morphological classifi¬ 
cations, by assuming there is no error in kinematic classi¬ 
fications. Only 80% (179/224) of the FRs were classified 
SO, so the probability of mistakenly classifying a “true” 
lenticular (i.e. an FR) as an E is p(E\FR) = 0.2. Of all 
the SRs, 65% (23/36) were classified as Es, so the prob¬ 
ability of mistakenly classifying a “true” elliptical (i.e. an 
SR) as an SO is p(S0\SR) = 0.35. Assuming ETGs are ei¬ 
ther Es or SOs (no misclassification of SRs or FRs as LTGs 
and vice versa), we can infer p(50|Ei?) and p(E\SR), as 
p(E\FR) +p(S0\FR) = p{E\SR) + p(S0\SR) = 1. Using 
these probabilities we can then estimate the number of mor¬ 
phological Es and SOs given the true number of FRs and 


SRs, 


where 


n(S0) 

n(E) 


n(FR) 

n(SR) 


p(S0|FR) p(S0|SR) 
p(E|FR) p(E|SR) 


(1) 


( 2 ) 


We can also invert Eq. □to estimate the number of FRs and 
SRs from the number of morphological Es and SOs, 


' n(FR) ' 

— p -1 

' ra(S0) ' 

n(SR) 


. "(E) . 


Using these equations, we confirm that the apparent 
(E:ETG) ratio would be 27% if the (SR:ETG) fraction is 
only 15% (and visa versa). Accordingly, we can appreci¬ 
ate why recent studies of the kT — E have found that the 
(SR:ETG) fraction remains roughly constant at 15%; in¬ 
deed, it could have been anticipated if the independence of 
the (E:ETG) ratio had been known a priori. 


4.3 The highest E-fractions are found in 
lower-density clusters 

We can use Eq. [3] to estimate the true number of FRs and 
SRs in individual clusters. Of great interest are the clus¬ 
ters with high (E:ETG) ratios, which are only found at low 
log(E 3 >. In Fig. [6] we see that there are 7 clusters with 
(E:ETG) >40% and log(Es)<1.5 (there are 26 clusters with 
(E:ETG) <0.4 and log(E 3 )<1.5). Such high (E:ETG) ratios 
are a result of much higher E-fractions and marginally lower 
SO-fractions: mean E-fractions are 25% and 12% for high and 
low (E:ETG) clusters, respectively; mean SO-fractions are 
32% and 45% for the high and low (E:ETG) clusters, respec¬ 
tively. The mean S+I fractions are very similar for both: 40% 
and 42% for the high and low (E:ETG) clusters, respectively. 
Consequently, it appears that the high (E:ETG) ratios are 
predominantly driven by a doubling of the E-fractions. Note 
that (E:ETG) ratio of 40% suggest (SR:ETG) ratio of 45% 
using Eq. [3] similarly an (E:ETG) ratio of 50% suggests an 
(SR:ETG) ratio of almost 70%. 

A key question to address is whether we have selected 
the upper envelope of a distribution (with large intrinsic 
scatter or large measurement uncertainty), or if these high 
(E:ETG) clusters represent a genuinely different type of clus¬ 
ter, perhaps at a different stage of cluster evolution (e.g. un¬ 
relaxed or merging). The scale of the random uncertainties 
above and below log(S 3 ) = 1.5 are comparable, so the high 
scatter at lower log(S 3 ) does not originate from increased 
measurement uncertainty In Fig. [6j there is an apparent gap 
between low (E:ETG) clusters and high (E:ETG) clusters 
below log(Ea) = 1.5. With the present data, we are unable to 
determine if this represents a genuine split and bimodality in 
the (E:ETG) ratio, or whether the distribution is continuous 
but sparsely sampled. 

In Fig. □ we showed that the high (E:ETG) clusters 
show a normal f(SI)-E relation, while the f(E)-E and f(S0)- 
E relations are offset from the usual trends and almost con¬ 
verge. This is not trivial: several studies of the kT — E re¬ 
lation have shown that in unrelaxed clusters, the SRs are 
not always found in the densest regions (Scott et al.||2014 
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. Furthermore, if the most massive galax- 
Es and/or SRs) merge in cluster merg¬ 
ers, then one might expect the high (E:ETG) clusters to be 
unrelaxed and show irregular T — E relations; this is not 
the case, they do show f(E)-E and f(S0)-E relations with 
higher numbers of ETGs in the denser environments. We 
therefore conclude that there is an intrinsic variation in the 
morphological content of lower-density clusters: some have 
double the normal E fractions and presumably much higher 
SR fractions. 


when forming progressively more massive clusters from more 
massive sub-units, there may be a process which dissipates 
the orbital angular momentum to increase the density of the 
galaxies, such as dynamical friction. Equally, as the most 
massive over-densities are the first to collapse, higher local 
densities in massive clusters may reflect the higher densities 
of galaxies at the formation epoch; lower mass clusters may 
only have virialised recently, after considerable expansion of 
the Universe. Both these possibilities rely on a link between 
global density and mass. 


Fogarty et al.J2014 
ies (predominantly 


4.4 The significance of the log(S 3 }— logE CT relation 


We have demonstrated that the E-fraction within clusters 
depends on the average projected local density of galaxies 
in those clusters, while the average E-fraction in the ETG 
population is independent of it. However, it is not imme¬ 
diately obvious what our ‘average projected local density’ 
measures. To address this problem, we investigated links be¬ 
tween log(Es} and a global measure of galaxy density. Fig. [8] 
illustrates that clusters with a higher average projected local 
density also have a higher projected global density. This is 
not immediately obvious, but equally perhaps not surpris¬ 
ing. What is peculiar is that when we define low-, medium- 
and high-density clusters based on this global measure of 
density, we do not recover the different f(E)-S relations seen 
in the lower panel of Fig. [4] Clearly, there is too much uncer¬ 
tainty in a global measure of density to recover these subtle 
changes. 

Our method of determining a global density relies on 
defining a single circular aperture to encompass the cluster 
galaxies, but many clusters are not circularly symmetric on 
the sky. This is an important limitation of the global ap¬ 
proach to measuring density. In Fig. [8j at a given average 
local density, there is a tendency for the global density to be 
lower when the circular aperture contained a smaller frac¬ 
tion of galaxies (i.e. the cluster was not well represented by 
the single centroid and circular aperture). This appears to 
be the dominant source of uncertainty (intrinsic scatter) in 
the log(E 3 )-logE CT relation, and may explain why we don’t 
recover the different f(E)-E relations when using global den¬ 
sity. 


Clusters frequently show multiple sub-structures and 
are not all circularly symmetric; for such ‘irregular’ mor¬ 
phologies, the original use of local projected density in D80, 
and the use of average local density here, have lead to a bet¬ 
ter understanding of how morphology links to environment. 
While the log(E 3 }- logE CT relation shows that there is a link 
between local environment and global environment, the sig¬ 
nificant scatter warns that we are not yet able to measure 
global environment reliably for all cluster morphologies. 

The log(E 3 ) logE CT relation shows some interesting fea¬ 
tures in its own right. There is a tendency for denser clusters 
to have a larger fraction of galaxies in the adaptive aperture, 
suggesting they are either more circularly symmetric and/or 
more concentrated in their galaxy distributions. The fact 
that there is a log(E 3 )-logE CT relation - a link between lo¬ 
cal density and global density - is not immediately obvious. 
For comparison, recall that more massive galaxies do not 


necessarily have higher stellar densities (see Fig. 1 of Cap 
pellari et alT|[2013j). It is not clear what mechanism could 
be responsible for this relation: in a hierarchical framework, 


4.5 Implications for formation mechanisms of Es 
and SOs 


The transformation of spirals into lenticulars is a common 
explanation for the fall/rise in the S+I/SO fractions with lo¬ 
cal density. However, the details of such a mechanism are far 


from certain and there are a number of contenders (Boselli 


& Gavazzi 20061. Recent evidence suggests that the ma¬ 


jority of cluster SOs were “pre-processed” and quenched in 
the group stages, before entering a cluster environment (e.g. 
|Mihos|[2004) Moran et al. 2007), via galaxy-galaxy interac¬ 


tions such as tidal interactions, harassment (Moore et al. 
19961 or gravitational heating (Khochfar fe Qstriker||2008 1. 


Indeed, the kT—T, relation of |Cappellari et al. | ( | 201 2[ ) probes 
mainly field and group environments: as they highlight, find¬ 
ing such a tight correlation in these lower-density environ¬ 
ments suggests that the T — E relation is initiated in the 


group stages. Similarly, Helsdon & Ponman (2003) find that 


the morphological segregation in X-ray bright groups is simi¬ 
lar to that of clusters, further supporting the suggestion that 
the group environment is mostly responsible for the T — E 


relations. Dressier et al. (2013) provide further support for 


pre-processing: they studied the stellar populations and sub¬ 
structures in intermediate redshift clusters and found a cor¬ 
relation between the passive and post-starburst fraction and 
the mass of the substructures, indicating the substructure 
environment (and not the overall cluster environment) dic¬ 
tates the quenching efficiency. Galaxy-galaxy interactions in 
the group stage are also far more likely to be related to local 
number density, which may explain its power in constructing 
the T — E relations. Cluster-scale effects like ram-pressure 
stripping ( Gunn fc Gott||l972| ), although evident ( |Chung| 
et al.|2009l, are not commonly thought to be the dominant 


route for creating the cluster SO population and the T— E re¬ 
lation (indeed, this was also the original conclusion of D80). 

So where do Es fit in to this picture? What distin¬ 
guishes Es from SOs morphologically, is the complete ab¬ 
sence of a disk. One can imagine numerous scenarios where 


a disk is destroyed. Historically, major mergers (Holmberg 


1941 Toomre & Toomre 

1972) were commonly used to ex- 

plain the assembly of Es l 

Zwicky|1959: Barnes|1988| Barnes[ 


& Hernquist 19926. But the major-merger scenario is not 


without its problems: unless collisions are head-on or the 
internal angular momentum of the stars aligned to cancel 
the orbital angular momentum of the merger, the remnants 
rotate too quickly and ar e too flattened to mimic slowly r o- 
tating elliptical galaxies ( White|l979a|b Bois et al. 1 2011 


4 Note that the simulations of|Bois et al. 

(]201l|) included a dark 

matter halo, while those of White (1979a 

b) did not. 
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A (multiple) minor merger scenario is now more commonly 
invoked to generate the low angular momentum state of the 


slowly rotating ellipticals (Khochfar et al. 2011b I. While this 
scenario has an appealing synergy in explaining the apparent 
size evolution of ETGs (e.g. Trujillo et al.|2007 Naab et al. 
2009), such a protracted assembly over time from many low- 


mass systems seems difficult to square with the high metal- 
licity and alpha-enhancement of the stellar populations in 
massive galaxies (implying high redshift star formation over 
short timescales, |Thomas et al.| |2005^ not to mention the 
abundance of Es at high redshift (Dressier et al.|1997 Post- 


man et aL||2005 Smith et al.|2005|. 

We have shown that Es obey a fundamentally differ¬ 
ent relation with local density compared to the disk sys¬ 
tems (SOs and S+Is): they reside in regions of the highest 
relative density. They either formed there, or accumulated 
there via some process. Take for example the Coma clus¬ 
ter: two super-giant, massive Es inhabit the centre of this 
cluster and they will merge in several Gyrs ([Gerhard et al. 
2007). These galaxies didn’t form in their present locations: 


dynamical friction from the cluster-scale dark matter halo 
has played a role in dissipating angular momentum to bring 
these two super-giant Es together. Significant sub-structure 
(i.e. galaxies) is associated with each of these super-giant 
Es and as in the general 3-body problem, the lower mass 
sub-structures may help transfer angular momentum away 
from the two most massive galaxies by being ejected. This 
will advance the merger, as may tidal effects. But the ejec¬ 
tion of substructures during the merger process will result 
in less of the initial orbital angular becoming locked into 
the remnant, and a lower stellar angular momentum state 
(a slow rotator). In addition, the substructure ejected on 
radial orbits could return, being stripped by tidal effects 
and even resulting in minor mergers to further lower the 
net angular momentum. But for the vast majority of clus¬ 
ter galaxies, mergers are rare (note merger signatures are 
long-lived, Sheen et al.||2012| Yi et al.]|2013 1. However, ha¬ 
rassment or high-velocity encounters ( jMoore et al.|1996 ) are 
common. All that is required to transform an SO to a E is 
destruction of the visible disk. Tidal effects and harassment 
may have the power to do this. The fall in the SO-fraction 
at the highest densities in Fig. [4] coincides with the peaks in 
the E-fractions, suggestion SOs are being transformed into 


Es there. Furthermore, D’Eugenio et al. (20151 have found 


that the ellipticities of passive galaxies decrease towards the 
centres of clusters, suggesting that the disks may be being 
disrupted or thickened as they fall into the denser regions 
of the clusters. Thus while the central SRs may form in the 
group stages, the remaining Es may have ‘formed’ (assumed 
their current morphology, previously being SOs) in the clus¬ 
ter itself. 

In [ Dressier et ah] (1 997[ ), studies of intermediate red¬ 
shift (z ~ 0.5) clusters yielded similar or higher E-fractions, 
with a deficit of SOs and an abundance of Ss. Indeed, typ¬ 
ical (E:ETG) ratios at z ~ 0.5 were as high as 75%, far 
higher than even the most E-rich clusters at low-z. Similar 
results were found at z~l, with an increased ETG-fraction in 
denser environments but similar ETG-fractions in the held 


(Smith et al. 2005 Postman et al. 2005 I. The conclusion 


from these observations is that the formation epoch of clus¬ 
ter Es is much earlier than cluster SOs, which formed later 
from quenched spirals. Accordingly, the clusters we see at 


low-z showing the highest (E:ETG) ratios in Fig. [6] may be 
relics of the earlier Universe which have undergone slower 
evolution to the present day. However, the implicit assump¬ 
tion that Es cannot be transformed into SOs (and visa versa) 
conflicts with the formation histories in |Naab et al. (20141, 
where many SRs (admittedly only 64% Es) are ‘formed’ (i.e. 
transform into their current kinematic morphology) at z< 1. 
Given that transformations between kinematic morphologies 
are seen in both directions in the simulations (FR «-» SR), 
transformations between visual morphologies may also be 
possible. 

Kinematic and visual morphology-density relations for 
intermediate redshift clusters may help us better understand 
this problem. The kinematic imprint of a disk is longer lived 
than the visible disk itself: comparing the SO and FR frac¬ 
tions at intermediate and low redshift will reveal if SOs are 
being transformed into (fast rotating) Es: a fall in the SO 
fraction and a rise in the E-fraction at higher densities for 
just the FR population would be interesting indeed. There 
is a school-of-thought that the high-redshift Universe was 
completely different and substantially more chaotic to the 
one we see locally. With mergers orders of magnitude more 
frequent (Fakhouri et al. 20101, star formation greatly en¬ 
hanced (Madau et al. 1996) and AGN activity similarly 
boosted (Cowie et al. 20031, we might well expect to create 
more dispersion-dominated structures with no stable disk. 
The same combined study of the kT — E and T — E relations 
could determine whether Es (or SRs) were really formed in 
their present structures at higher redshifts; or if the central 
galaxies in clusters gradually form their low angular momen¬ 
tum state from minor mergers. We can further determine if 
the statistics of our visual morphologies at higher redshift 
tally with those at low redshift, using the S0:FR and E:SR 
ratios, which is crucial given the doubts raised by | Andreon| 
(19981 and Holden et al. (20091. 


5 CONCLUSIONS 


We revisited the original data of Dressier (1980aJ, updat¬ 
ing the cluster redshifts and projected densities. We use the 
average local projected density to characterise clusters into 
low-, medium- and high-density classes and show that: 


• the rise in elliptical fraction with local density for low- 
density clusters occurs at lower densities than the same rise 
for high-density clusters; thus, the distribution of elliptical 
galaxies within a cluster depends on the relative local pro¬ 
jected density of galaxies, which is not the case for the lentic¬ 
ular or late-type fractions; 

• the (E:ETG) ratio in a cluster is independent of the av¬ 
erage local density of galaxies in that cluster, with a constant 
value of around 30%; in the fast and slow rotator paradigm, 
this is consistent with a constant (SR:ETG) fraction of 15%; 
conversely, the (S0:Disk) ratio varies strongly with average 
local density; 

• contrary to the canonical T — E relation, clusters with 
the highest (E:ETG) ratios (and overall elliptical fractions) 
are found in lower-density clusters; the T — E relations for 
these high (E:ETG) ratio clusters show a convergence of the 
f(E)-E and f(S0)-E relations; 

• there is a simple relation between the average local pro¬ 
jected density of galaxies and the global density of galaxies 
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in a cluster, such that clusters with a higher average lo¬ 
cal density have a higher global density; however our global 
measure of density remains crude and poorly characterises 
irregular cluster morphologies. 
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are associated with the z«0.03 over-density compared to the 
0.45<z<0.75 structures. 

This paper has been typeset from a TgX/ DTjgK file prepared 
by the author. 


APPENDIX A: PROJECTION EFFECTS IN 
ABELL 838, ABELL 1631 & DC2349-28 


While updating the cluster redshifts using the latest values 
from NEEQ we further verified the cluster redshifts using 
the galaxy redshifts. Galaxy redshifts were compiled from: 
NED (mostly from the 6dF survey, Jones et al.||2004 l; the 
2MASS redshift survey (Huchra et al. 20121, the WINGS 
survey (|Moretti et al. 1 12014| ) ; the CAIRNS survey (Rines 


et al.|2003|). We found that the NED redshifts for the clus¬ 


ters Abell 838, Abell 1631 and DC2349-28 (0.0502, 0.01394 
& 0.0648, respectively) were inconsistent with the redshifts 
of all the galaxies. A detailed investigation of the redshift 
distributions for each revealed that they are each most likely 
projections of two or more clusters along the line of sight. 

For Abell 838, only 19 of the 62 galaxies were success¬ 
fully cross-matched with the aforementioned redshift sur¬ 
veys. For these 19 galaxies, the redshift distribution is bi- 
modal and equally distributed at z«0.023 & z«0.052. How¬ 
ever, when considering all galaxies in the spectroscopic cat¬ 
alogues within a radius of 0.8°, more than twice as many 
galaxies are associated with the z«0.023 over-density than 
the one at z«0.052. As such, we associate Abell 838 with 
a cluster at z=0.023. For Abell 1631, we successfully cross¬ 
match 124 of the 139 galaxies in D80 and find that only 17 
are associated with an over-density at z=0.014, while 103 are 
associated with an over-density at z«0.047. Accordingly, we 
associate Abell 1631 with a cluster at z=0.047. Finally, for 
DC2349-28, we cross-match 55 of 68 galaxies and find 28 
galaxies associated with an over-density at z~0.03, and 22 
galaxies associated with a number of smaller over-densities 
at 0.045<z<0.075. When considering all galaxies with red¬ 
shifts within a radius of 0.93 °, we find that twice as many 


5 http://ned.ipac.caltech.edu 


















